Theory of swimming filaments in viscoelastic media 
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Motivated by the swimming of sperm in the non-Newtonian fluids of the female mammalian repro- 
ductive tract, we examine the swimming of fllaments in the nonlinear viscoelastic Upper Convected 
Maxwell model. We obtain the swimming velocity and hydrodynamic force exerted on an infinitely 
long cylinder with prescribed beating pattern. We use these results to examine the swimming of a 
simplified sliding-filament model for a sperm fiagellum. Viscoelasticity tends to decrease swimming 
speed, and changes in the beating patterns due to viscoelasticity can reverse swimming direction. 
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The physical environment of the cell places severe con- 
straints on mechanisms for motility. For example, viscous 
effects dominate inertial effects in water at the scale of a 
few microns. Therefore, swimming cells use viscous resis- 
tance to move, since mechanisms that rely on imparting 
momentum to the surrounding fluid, such as waving a 
rigid oar, do not work [H, [3]. The fundamental princi- 
ples of swimming in the low-Reynolds number regime of 
small-scale, slow flows have been established for many 
years 0, H, 0, 0|i yet continue to be an area of active 
research. However, when a sperm cell moves through the 
viscoelastic mucus of the female mammalian reproductive 
tract, the theory of swimming in a purely viscous fluid is 
inapplicable. Observations of sperm show that they are 
strongly affected by differences between viscoelastic and 
viscous fluids. In particular, the shape of the flagellar 
beating pattern as well as swimming trajectories and ve- 
locities depend on the properties of the medium [1, 0, S] ■ 

The interplay of medium properties and flagellar motil- 
ity or transport also arises in other situations, such as 
spirochetes swimming in a gel Q, and cilia beating in 
mucus to clear foreign particles in the human airway [lo| . 
Motivated by these phenomena, we develop a theory for 
swimming filaments in a viscoelastic medium. We begin 
by analyzing the swimming of an infinite filament with 
a prescribed beating pattern in a fluid described by the 
Upper Convected Maxwell (UCM) model [uj. We de- 
duce the hydrodynamic force per unit length acting on 
the filament and the swimming velocity to leading order 
in the deformation of the filament. Our results extend 
the findings of Lauga 0, who considered a variety of 
fading memory models for the case of a prescribed beat 
pattern on a planar sheet. We further apply our results 
to a simple model fiagellum with active internal forces, 
and find that changes in fiagellum shapes play a crucial 
role in distinguishing the effects of viscoelastic media. 

Newtonian fluids are characterized by a simple consti- 
tutive relation, in which stress is proportional to strain 
rate. Non-newtonian fiuids cannot be characterized by a 
simple universal constitutive relation, and exhibit a range 
of phenomena such as elasticity, shear thinning, and yield 



stress behavior. We choose to focus our attention on flu- 
ids with fading memory, in which the stress relaxes over 
time to the viscous stress. We consider small amplitude 
deflections of an infinite filament of radius a (Fig. [1^). 
Since the swimming velocity of a filament is second order 
in the deflection amplitude 0, 0] , linear models for fad- 
ing memory such as the Maxwell model are insufficient 
for studying swimming (l^ . 

Therefore, we use the simplest nonlinear constitutive 
relation incorporating elastic effects, the UCM model. 
This model is appropriate for a polymer solution in which 
the viscosity of the Newtonian solvent is disregarded: 



T + Xt — 777. 



(1) 



In this equation r is the deviatoric stress, A is the relax- 
ation time, T = dtT + v- Vr — (Vv)^-t — r- Vv is the 
upper-convected time derivative of r, v is the velocity, rj 
is the polymer viscosity, and 7 = Vv-l-(Vv)^ is the strain 
rate. The nonlinear terms of the upper-convected deriva- 
tive make the constitutive relation insensitive to trans- 
lational and rotational motion of material elements \11 \ . 



The UCM fluid responds as an elastic solid when sub- 
ject to a rapidly varying stress, and as a viscous liquid 
when subject to a slowly varying stress. When A = 0, 
the constitutive relation ([T]) is Newtonian. Since inertia 
is unimportant, the motion of the medium is governed 
by force balance, — Vp + W ■ t — 0, or 



(l + A9t)Vp 



AV-T, 



(2) 



where T = v-Vt — (Vv)^t — t-Vv, and we have assumed 
incompressibility, V • v = 0. 

To calculate the swimming velocity of the filament, we 
prescribe a beating pattern that is independent of load, 
and solve ^ for the fiow, imposing no-slip boundary 
conditions at the surface of the filament. Material points 
on the the surface of the filament are parameterized by 
z and ip as in Fig. [1^. The positions of points on the 
surface are given by 



r(z, ip, t) = [h{z, t) + a cos{tp)]ic + a sin(-0)y - 



(3) 
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FIG. 1: a) Cylinder with imposed traveling wave of trans- 
verse displacements in a viscoelastic fluid, b) Active flagel- 
lum model. The active elements inside the flagellum represent 
motors that slide the filaments relative to each other. 



where h{z, t) — Re^^ ^ hq^ exp{iqz — iujt)). Note that 
we work in the frame in which material points of the 
cylinder move in planes of constant z. We will find that 
the no-slip boundary conditions can only be satisfied if 
there is a uniform flow along the z axis. In the lab frame 
in which the fluid is at rest at infinity, the flow corre- 
sponds to the swimming velocity of the cylinder. 

Invoking small amplitude, long wavelength distortions 
for simplicity, we solve the equations order by order in 
the displacement amplitudes qhq^^ to lowest order in 
f/logga, and assuming that hq^^/a ^ f 
the dynamic variables as expansions in qhqi^, 



Expressing 
such as 

V — v'^' -|- v^^-* -t- the no-slip boundary conditions 
in cylindrical coordinates (r, 0, z) up to second order are 



Ax = v(i) 



+ v^^^ + h (cos ijj r - sin ijj (i>j ■ Vv^^^ (4) 



where h is evaluated at and v, f, and at 

{a,tjjjZ,t). We used the fact that the azimuthal angle 
(j) of the material point labeled ip on the cylinder surface 
is given hy (j) ~ ip — (h/a) sinip, and the radius of this 
point is given by r w a -I- hcosip. There is no first order 
contribution to T and the first order dynamical equations 
are 



(f + Aai)r« = 777 



.(1) 



{l + \dt)Vp 



(1) - 



2^(1) 



(5) 
(6) 



with V • v^^) = 0. Since Eq. ^ is Stokes equation with 
a modified pressure, the solution is readily found ^IS^, 
and the first order flow v*^^) is the same as in the purely 
viscous case, in the limit qa <^1. The formulas for this 
flow are given by Taylor in |4j. Just as in Taylor's case, 
the swimming velocity vanishes to first order in qhq^. 
Using the pressure and stress, wc find that the force per 
unit length is purely in the x direction: 

C(z,t) = Re 5] 



q,uj 



(1 — iXuj) logga 



(7) 



This expression is consistent with the results of Fulford, 
Katz, and Powell, who used resistive force theory for a 
filament beating in a linear Maxwell fluid to show that 
there is no change in the swimming velocity relative to 



the viscous case, even to second order in qhq^^ 

In our problem, however, the nonlinearities make the 
viscoelastic swimming velocity different from the viscous 
swimming velocity. To second order, 

(1 + A9t)T(2) = ^7(2) _ at(2) (8) 
(1 -f A9t)Vp(') = r;V2v(2) - AV-T(2), (9) 

where T*^2) be calculated using only the flrst order 
stresses and velocity flelds. To flnd the time-averaged 
swimming velocity, we need only examine the velocity 
fields averaged over time and 0. To second order, we find 
that there is a uniform flow at infinity, corresponding to 
a swimming velocity in the lab frame of 



U(2) 



1 - 

1 

q,uj 



\hq^\'^quj 



(10) 



For a single traveling wave the direction of swimming 
is opposite the direction of motion of the traveling 
wave. Our result, valid to first order in an expansion in 
1/ log(qa), is precisely the same as the case of a traveling 
wave on a planar sheet [13] ■ 

We have dealt with the effects of viscoelasticity on a 
swimming filament with prescribed shape h{z, t). In con- 
trast, the beating patterns and the swimming velocity of 
real sperm are affected by the medium, suggesting that 
prescribing the shape changes of the swimmer may miss 
important effects. Therefore, we consider a sliding fila- 
ment model in which we prescribe active internal bending 
forces 15|, |l6| , and solve for the flagellum shape as well 
as the swimming speed. Our task is greatly simplified 
since the shape is required to an accuracy of first order 
in deflection only. The model is shown in Fig. [Us. We 
continue to assume the flagellum has a cylindrical cross- 
section as in Fig[T^ and deflection from the z-axis /ix, but 
now we assume the flagellum has finite length L and two 
incxtensible longitudinal filaments with constant lateral 
spacing 2a. Motors along the flagellum attach to both 
filaments and tend to slide them past each other. Sliding 
is prohibited at the end near the head, which is omitted 
for simplicity. Assuming a planar shape of the flagellum, 
the moment M = My acting at a cross-section of the 
flagellum consists of a passive resistance to bending, and 
an active part due to the sliding motors. To flrst order 
in deflection /i, 



M{z) = Ah" -2a f^dz, 



(11) 



where A is the bending stiffness, primes denote deriva- 
tives with respect to 2, and /m is the force per unit length 
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that the lower filament of Fig. [Jd exerts on the upper fil- 
ament. For simplicity we disregard any elastic or viscous 
effects arising from proteins linking the filaments. 

The balance of internal forces and hydrodynamic forces 
determines the instantaneous shape of the flagellum. To 
find the internal force per unit length /i„t, consider mo- 
ment balance on an element of the flagellum to first or- 
der in deflection [l3|: M' + N — Q, where N is the shear 
force acting on a cross section in the x direction. Using 
/int = N' yields /int = -Ah"" - 2af^,. We choose a 
sliding force fm = Re[/exp(zA:z — iojt)]. To linear order, 
the shape change occurs with only one frequency w, and 
h{z,t) = Re[h{z) exp{—iiut)]. 

Since the internal forces /int are expressed in real space, 
it is convenient to write the hydrodynamic force Eq. ([7|) 
in real space. The logarithm in Eq. ([7|) varies slowly with 
q and we replace it with a constant, log(ga) ~ log(a/L), 
as is commonly done in resistive force theory. Thus 

ffluid - Re^-^ (-zc^/i(z)e-*-*) X, (12) 

where C± ~ 4:nrj / log {L / a) . 

We non-dimensionalize our equations of motion by 
measuring lengths in terms of L, in terms of 
A/(2aL^), and time in terms of w^^. For notational sim- 
plicity, after scaling we use the same symbols for the new 
quantities. The equation of motion for the active flagel- 
lum is fint + ffluid = 0, or in non-dimensional form, 

~'^P^ h + h"" - ikf ^ 0. (13) 
1 — iDe 

The dimensionless "Sperm number" Sp = L{lij(±/AY^^ 
is the fourth root of the ratio of the bending relax- 
ation time of the flagellum to the period of the trav- 
eling wave, and the Deborah number De = Xto mea- 
sures the importance of elastic effects. Following [l6| . 
we estimate L = 40 ^m, a = 20 nm, A = A x 10"^^ N- 
m^, tj/(27r) = 30s-i, and C_l = 2 x 10"^ N-s-ni^^ 
and a dimensional magnitude of internal sliding forces 
/m = 4pN/(24nm). Therefore, Sp « 7, and the dimen- 
sionless amplitude / w 13. For sperm which oscillate at 
frequency 25-50 Hz in cervical mucus with time constant 
A « 1-10 s [11], we take De « 100. Equation ^ must 
be supplemented by boundary conditions. For simplicity 
we forbid transverse motion of the head h{0) = 0, and 
suppose the connection between the head and the flagel- 
lum cannot support a moment: h"{0) + /in(z)dz = 
(see [l^ for a more realistic treatment of the moment 
boundary condition). The boundary condition at the 
other end is zero force, —h"'{l) + /m(l) = 0, and zero 
moment, h"{l) — 0. 

Representative plots of the beating patterns are shown 
in Fig. [21 Two lengthscales are apparent. The first is 
the wavelength of the sliding force, 27r/fc, readily ap- 
parent in the change of beating patterns from fc = 



to fc = Stt/L at Sp = 7. The second lengthscale ^ arises 
from the interaction of bending and hydrodynamic forces: 
C/L = |1 -iDc|i/4/Sp. The flagellum behaves like a rigid 
rod for large £,/L, and is floppy for small £,/L. Changing 
the viscoelastic properties of the fluid affects the beating 
shapes through ^, as can be seen by comparing beating 
patterns for De = and De = 100 in the flagellum with 
k = Stt/L and Sp = 7. In Fig. [51 the amplitude of the 
sliding force is selected so that the maximum displace- 
ment of the flagellum is L/10. Since ^ increases with De, 
smaller driving forces are required to produce the same 
amplitude motion. Besides affecting ^, viscoelasticity af- 
fects phase behavior: for De = 0, we obtain traveling 
wave beating patterns, while for large De, in the elastic 
limit, we obtain standing wave beating patterns. 

By inserting the beating pattern hq^ into Eq. (jlOp . we 
can calculate the swimming speed of a flagellum with 
prescribed internal forces. Although our flagellum is fl- 
nite, we may apply Eq. pUj) to our calculated beating 
pattern by periodically replicating it to infinite extent in 
z. Physically, we are ignoring end effects. It is useful to 
rewrite Eq. (|10p in real space, for an infinite flagellum 
with period L and single oscillation frequency uj: 

IL[1 + UC j Jo 

where the dot denotes the time derivative and the brack- 
ets denote time-averaging. The swimming velocity is in- 
dependent of whether we consider the period L or nL; 
thus, the periodically replicated flagellum has the same 
velocity as the single flagellum, up to end effects. 

The velocities of active flagella with sliding forces of 
varying wavelength are plotted in Fig. [31 Viscoelas- 
tic effects can produce qualitative changes in the swim- 
ming velocity. Most notably, as De is increased, the ve- 
locity can change direction, for example in the case of 
k = Stt/L. When k = Stt/L and A = 0, there is a 
small amplitude traveling wave moving in the positive 
z-direction with wavelength w L/4 driven by the slid- 
ing forces, but there is also a larger amplitude traveling 
wave moving in the negative direction (Fig. [5]). A fila- 
ment swims in the direction opposite to traveling waves 
moving in the filament so these push the flagellum in op- 
posite directions, but the net velocity is in the positive 
direction. As De increases, the larger amplitude wave be- 
comes a standing wave, and only the smaller amplitude 
traveling waves remain, causing the fiagellum to move in 
the negative direction. We emphasize that this effect is 
not due to the (l-fDe^) correction factor to the swimming 
speed in Eq. (|14p , which can never reverse the swimming 
direction. For swimming in viscoelastic fluids it is crucial 
to allow the flagellum shape to respond to changed forces 
rather than prescribe fixed beating patterns. 

We conclude by indicating directions for future work. 
More realistic constitutive relations should be explored. 
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FIG. 2: (color). Dimensionless amplitude h/L versus dimensionless arclength s/L for beating flagella with internal sliding 
forces. A half-cycle (red, orange, green, light blue, blue) of the pattern is shown for viscous (De = 0) and viscoelastic (De = 100) 
cases. For Sp = 7, we show beating patterns for internal sliding forces with k varying from (uniform force) to Stt. We also 
show how changing the sperm number affects the shape for k — Sn/L. On the plots for = 0, we show ^. At the top of each 
plot, the (dimensionless) magnitude / required to produce motion with amplitude O.IL is shown. 
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FIG. 3: Dimensionless swimming velocities versus De, for 
Sp = 7 and maximum flagellum displacement O.IL, for fc = 
(long dashes), k — 3n/L (short dashes), k = Svr/L (solid). 
The inset shows detail for the same plot demonstrating that 
viscoelastic effects can reverse the direction of swimming. 



since the UCM model is invalid for sufficiently high ex- 
tensional flows and does not include shear thinning ef- 
fects. Realistic modeling should also include the effects of 
ends and large displacements disregarded here. Finally, 
it is important to determine whether the forces imposed 
on the filaments by motors are the same in viscoelastic 
and viscous media, or, instead, sperm respond to changes 
in load. We note that the observation that beating fre- 
quency is dependent upon the medium 0, S| suggests 
the latter. 
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